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Before the launch of the Fermi satellite only two classes of Active Galactic Nuclei (AGN) were known to generate 
relativistic jets and thus to emit up to the 7-ray energy range: blazars and radio galaxies, both hosted in giant 
elliptical galaxies. The first four years of observations by the Large Area Telescope (LAT) on board Fermi 
confirmed that these two populations represent the most numerous identified sources in the extragalactic 7-ray 
sky, but the discovery of variable 7-ray emission from 5 radio-loud Narrow-Line Seyfert 1 (NLSyl) galaxies 
revealed the presence of a possible emerging third class of AGN with relativistic jets. Considering that NLSyl 
are thought to be hosted in spiral galaxies, this finding poses intriguing questions about the nature of these 
objects, the knowledge of the development of relativistic jets, and the evolution of radio-loud AGN. In this 
context, the study of the radio-loud NLSyl from radio to 7-rays has received increasing attention. Here we 
discuss the radio-to-7-rays properties of the 7-ray emitting NLSyl, also in comparison with the blazar scenario. 


1. Introduction 

Relativistic jets are the most extreme expression of 
the power that can be generated by a superluminal 
black hole (SMBH) in the center of an active galac¬ 
tic nucleus (AGN), with a large fraction of the power 
emitted in 7-rays. Before the launch of the Fermi 
satellite only two classes of AGN were known to gen¬ 
erate these structures and thus to emit up to the 7- 
ray band: blazars and radio galaxies, both hosted in 
giant elliptical galaxies 0- The first 4 years of obser¬ 
vation by the Large Area Telescope (LAT) on board 
Fermi confirmed that the extragalactic 7-ray sky is 
dominated by radio-loud AGN, being mostly blazars 
and some radio galaxies Q. However, the discovery 
by Fermi -LAT of variable 7-ray emission from a few 
radio-loud narrow-line Seyfert 1 (NLSyl) galaxies re¬ 
vealed the presence of a possible third class of AGN 
with relativistic jets [lj. 

NLSyl are a class of AGN identified by 143) and 
characterized by their optical properties: narrow per¬ 
mitted lines (FWHM (H/3) < 2000 kms" 1 ), [OIII]/H/3 
< 3, and a bump due to Fe II (e.g., 01). They also 
exhibit strong X-ray variability, steep X-ray spectra, 
relatively high luminosity, and substantial soft X-ray 
excess (e.g., (0). These characteristics point to sys¬ 
tems with smaller masses of the central black hole 
(10 6 TO 8 Mq) than blazars and radio galaxies, and 
higher accretion rates (close to or above the Edding¬ 
ton limit). NLSyl are generally radio-quiet (radio¬ 
loudness R < 10), with only a small fraction (< 7%; 
0]) classified as radio-loud. Objects with high val¬ 
ues of radio-loudness (R > 100) are even more sparse 
(~2.5%), while generally ~15% of quasars are radio- 
loud. Considering also that NLSyl are thought to 
be hosted in spiral galaxies, their detection in 7-rays 
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Figure 1: Fermi -LAT (E > 100 MeV) light curve of SBS 
0846+513 obtained during 2012 April 1 - August 28 with 
7-day (black circles) and 1-day (red triangles) time bins. 
Arrows refer to 2-cr upper limits. Upper limits are 
computed when TS <10. Adapted from [2ll ]. 

poses intriguing questions about the nature of these 
sources, the production of relativistic jets, the mech¬ 
anisms of high-energy emission, and the cosmological 
evolution of radio-loud AGN. 


2. The 7 -ray view of NLSyl 

Five radio-loud NLSyl galaxies have been de¬ 
tected at high significance by Fermz-LAT so far: 1H 
0323+342, SBS 0846+513, PMN J0948+0022, PKS 
1502+036, and PKS 2004-447 (J i, 0 , with a red- 
shift between 0.061 and 0.585. The average apparent 
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Figure 2: XMM-Newton EPIC pn (black), MOS1 (red) 
and MOS2 (green) data of PMN J0948+0022 shown as a 
ratio to a power law with T = 1.48. Adapted from [2li . 


while a soft X-ray excess is evident below ~2 keV 
(Fig. [2]). The origin of the soft X-ray excess is still an 
open issue both in radio-quiet and radio-loud AGN 
(e.g., (HI). Such a Seyfert component is a typical fea¬ 
ture in the X-ray spectra of radio-quiet NLSyl, but 
it is quite unusual in jet-dominated AGN, even if not 
unique (e.g., PKS 1510-089; (HI). In the case of PMN 
J0948+0022, the statistics did not allow us to dis¬ 
tinguish between different models for the soft X-ray 
emission. Models where the soft emission is partly 
produced by blurred reflection, or Comptonisation of 
the thermal disc emission, or simply a steep power- 
law, all provide good fits to the data. A multicolor 
thermal disc emission also gives a comparable fit, but 
the temperature is too high (kT = 0.18 keV) and is 
incompatible with a standard Shakura & Sunyaev ac¬ 
cretion disc HU . 


isotropic luminosity of these sources in the 0 . 1-100 
GeV energy band is between 10 44 erg s _1 and 10 4 ' 
erg s _1 , a range of values typical of blazars pj|. This 
may be an indication of a small viewing angle with 
respect to the jet axis and thus a high beaming factor 
for the 7 -ray emission, similarly to blazars. Several 
strong y-ray flares were observed from SBS 0846+513 
(Fig. [TJ) and PMN J0948+0022, reaching a peak ap¬ 
parent isotropic 7 -ray luminosity of ~ 10 48 erg s _1 , 
comparable to that of the bright FSRQ HI US H3| • 
In particular, SBS 0846+513 and PMN J0948+0022 
showed a 7 -ray flaring activity combined with a mod¬ 
erate spectral evolution [3 HU, a behaviour that 
was already observed in bright flat spectrum radio 
quasars (FSRQ) and low-synchrotron-peaked BL Lacs 
[2). Variability and spectral properties of these two 
NLSyls in 7 -rays indicate a blazar-like behaviour. An 
intense 7 -ray flaring activity was observed by LAT 
also from 1H 0323+342 [HI- This is another indica¬ 
tion that radio-loud NLSyl are able to host relativistic 
jets as powerful as those in blazars. 


3. X-ray properties 


The X-ray spectra of NLSyl are usually charac¬ 
terized by a steep photon index (r x > 2 , HU). 
On the contrary, a relatively hard X-ray spectrum 
was detected in the Swift/X RT observations of SBS 


0846+513 [3,120], PMN J0948+0022 [22 
1H 0323+342 [3, and PKS 1502+036 |17|. This sug¬ 
gests a significant contribution of inverse Compton ra¬ 
diation from a relativistic jet, similar to what is found 
for FSRQ. 

The high quality XMM-Newton observation of PMN 
J0948+0022 performed in 2011 May allowed us to 
study in detail its X-ray spectrum, as reported in [22]. 
The spectral modelling of the XMM-Newton data of 
PMN J0948+0022 shows that emission from the jet 
most likely dominates the spectrum above ~2 keV, 


4. Radio properties 

On pc scale a core-jet structure was observed for 
SBS 0846+513 [3 (Fig.®, PKS 2004— 447j4llL PKS 
1502+036 HU. and PMN J0948+0022 jUfS], al¬ 
though the jet in the two latter sources is significantly 
fainter than that observed in the former two sources. 
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Figure 3: 15.3 GHz VLBA image of SBS 0846+513. On 
the image we provide the peak flux density, in 
mjy/beam, and the first contour intensity (f.c., in 
mjy/beam) that corresponds to three times the noise 
measured on the image plane. Contour levels increase by 
a factor of 2. The beam is plotted on the bottom left 
corner of the image. Component W1 is the core region, 
W2 is a knot, and E is the jet structure. Adapted from 
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Figure 4: The separation between the core component of 
SBS 0846+513 and the knot ejected in 2009 as a function 
of time. The solid line represents the regression fit to the 
15 GHz VLBA MOJAVE data, while the dotted lines 
represent the uncertainties from the fit parameters. 
Dashed lines indicate the beginning and the peak of the 
radio flare. Adapted from [2ll |. 



Figure 5: SED and models for the 2013 and 2011 activity 
states of PMN J0948+0022. The filled circles are the 
data from the 2013 flaring state, and the open squares 
are the data from the 2011 intermediate state taken from 
[22} |. The dashed curve indicates the disc and coronal 
emission, and the dotted line indicates the thermal dust 
emission. Solid lines represent models consistent with 
scattering dust torus radiation, while the dashed-dotted 
curve represents a model consistent with the scattering of 
BLR radiation. Arrows refer to 2<r upper limits on the 
source flux. The VERITAS upper limits are corrected for 
EBL absorption using the model of [261. Adapted from 

El- 


The analysis of the 6 -epoch data set of SBS 0846+513 
collected by the MOJAVE programme during 2011- 
2013 indicates that a superluminal jet component is 
moving away from the core with an apparent angular 
velocity of (0.27+0.02) mas yr -1 (Fig.|4}, correspond¬ 
ing to (9.3±0.6)c |20j . This apparent superluminal 
velocity indicates the presence of boosting effects for 
the jet of SBS 0846+513. On the contrary, VLBA ob¬ 
servations did not detect apparent superluminal mo¬ 
tion at 15 GHz for PKS 1502+036 during 2002-2012, 
although the radio spectral variability and the one¬ 
sided jet-like structure seem to require the presence of 
boosting effects in a relativistic jet Hi- 

Strong radio variability was observed at 15 GHz 
during the monitoring of the OVRO 40-m telescope 
of PMN J0948+0022 plJp, PKS 1502+036 Hfj, 
and SBS 0846+513 [IflEcj. An inferred variability 
brightness temperature of 2.5xl0 13 K, 1.1 xlO 14 K, 
and 3.4xlO 11 K was obtained for PKS 1502+036, 
SBS 0846+513, and PMN J0948+0022, respectively. 
These values are larger than the brightness tempera¬ 
ture derived for the Compton catastrophe [Hi], sug¬ 
gesting that the radio emission of the jet is Doppler 
boosted. On the other hand, a high apparent bright¬ 
ness temperature of 10 13 K, comparable to that of the 
7 -ray NLSyl, was observed for TXS 1546+353. How¬ 
ever, no 7 -ray emission has been detected from this 
source, so far [423. Moreover, an intensive monitor¬ 
ing of these 7 -ray NLSyl from 2.6 GHz to 142 GHz 
with the Effelsberg 100-m and IRAM 30-m telescopes 


showed, in addition to an intensive variability, spec¬ 
tral evolution across the different bands following evo¬ 
lutionary paths explained by travelling shocks, typical 
characteristics seen in blazars jsj. 

5. Multifrequency variability and SED 
modelling 

The first spectral energy distributions (SED) col¬ 
lected for the four NLSyls detected in the first year of 
Fermi operation showed clear similarities with blazars: 
a double-humped shape with a first peak in the 
IR/optical band due to synchrotron emission, a sec¬ 
ond peak in the MeV/GeV band likely due to inverse 
Compton emission, and an additional component re¬ 
lated to the accretion disc in UV for three of the four- 
sources. The physical parameters of these NLSyl are 
blazar-likc, and the jet power is in the average range 
of blazars [ 1 ]. 

For PMN J0948+0022 we compared the broad¬ 
band SED of the 2013 flaring activity state with that 
from an intermediate activity state observed in 2011 
(Fig. [5]). Contrary to what was observed for some 
FSRQ (e.g., PKS*0537-441; (Hi) the SED of the 
two activity states, modelled as synchrotron emission 
and as an external Compton scattering of seed pho¬ 
tons from a dust torus, could not be modelled by 
changing only the electron distribution parameters. 
A higher magnetic field is needed for the high activity 
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Figure 6 : Multifrequency light curve for SBS 0846+513 
during 2011 December - 2013 January. The data sets 
were collected (from top to bottom) by Fermi -LAT 
( 7 -rays, 0.1 — 100 GeV; in units of 10“' ph cm -2 s _1 ), 
Swift-XRT (0.3—10 keV; in units of 10 -13 erg cm -2 s _1 ), 
Swift -UVOT (w2 and u bands; in units of /ijy) and 
OVRO (15 GHz; in units of mjy). Arrows refer to 30- 
upper limits on the source flux densities for the w 2 and u 
bands, and to 2 a- upper limits on the source fluxes for the 
7 -ray light curve. Adapted from 2(j- 


state, consistent with the modelling of different activ¬ 
ity states of PKS 0208—512 [l3| . We also modelled 
the 2013 flaring state assuming Compton-scattering 
of broad line region (BLR) line radiation. The model 
reproduces the data as well as the scattering of the 
IR torus photons. However, we note that the BLR 
scattering model requires magnetic fields which are 
far from equipartition. 

We also compared the SED of SBS 0846+513 during 
the flaring state in 2012 May with that of a quiescent 
state. Similar to PMN J0948+0022, the SED of the 
two different activity states, modelled by an exter¬ 
nal Compton component of seed photons from a dust 
torus, could be fitted by changing the electron distri¬ 
bution parameters as well as the magnetic field ( 20 ]. 
A significant shift of the synchrotron peak to higher 
frequencies was observed during the 2012 May flaring 
episode, similar to FSRQ (e.g., PKS 1510—089; (HI). 
Contrary to what is observed in PMN J0948+0022, 
no significant evidence of thermal emission from the 
accretion disc has been observed in SBS 0846+513. 

A complex connection between the radio and 7 -ray 
emission was observed for SBS 0846+513 and PMN 
J0948+0022, where 7 -ray and radio flares have not a 
similar behaviour, as discussed in detail in l 20 l l 22 l |28| . 


Optical intraday variability has been reported for 
PMN J0948+0022 by 32. 35] . [38| , sometimes asso¬ 
ciated with a significant increase of the optical polar¬ 
isation percentage, indicating a relativistic jet as the 
most likely origin for the optical emission in this ob¬ 
ject. 

At Very High Energy (VHE; E > 100 GeV), VER¬ 
ITAS observations of PMN J0948+0022 were carried 
out during 2013 January 6-13, after the 7 -ray flare 
observed by Fermi -LAT on 2013 January 1. These 
observations resulted in an upper limit of F>o.2TeV 
< 4xl0 -12 ph cm -2 s _1 [23j]. The lack of detection 
at VHE could be due to different reasons: 1) The dis¬ 
tance of the source (z = 0.5846) is relatively large and 
most of the GeV/TeV emission may be absorbed due 
to pair production from 7 -ray photons of the source 
and the infrared photons from the extragalactic back¬ 
ground light (EBL). However, we must note that the 
most distant FSRQ detected at VHE up to now, 3C 
279 Q is at a comparable distance. 2) The VERI¬ 
TAS observations were carried out a few days after 
the peak of the 7 -ray activity, thus covering only the 
last part of the MeV/GeV flare. 3) Considering the 
similarities with FSRQ, a BLR should be present in 
these NLSyl. The presence of a BLR could produce 
a spectral break due to pair production, suppressing 
the flux beyond a few GeV and preventing a VHE de¬ 
tection. However, the detection at VHE of the FSRQ 
3C 279, PKS 1510-089 @,0, and 4C +21.35 Q have 
shown that the spectrum of some FSRQ extends to 
VHE energies during some flares, indicating that the 
7 -rays may be produced outside the BLR during those 
high-activity periods. The same scenario may apply 
to the 7 -ray emitting NLSyl. 


6 . Radio-loudness, host galaxy, and jet 
formation 

The mechanism at work for producing a relativistic 
jet is not clear. In particular, the physical parame¬ 
ters that drive the jet formation is still under debate. 
By considering that NLSyl are thought to be hosted 
in spiral galaxies (e.g., |25|) the presence of a rela¬ 
tivistic jet in these sources seems to be in contrast to 
the paradigm that such structures could be produced 
only in elliptical galaxies (e.g., HI). The most power¬ 
ful jets are found in luminous elliptical galaxies with 
very massive central BH and low accretion rates (e.g., 
HHH). This was interpreted as an indirect evidence 
that a high spin is required for the jet production, 
since at least one major merger seems to be neces¬ 
sary to spin up the SMBH. At the same time, low 
accretion rates, which are associated with geometri¬ 
cally thick advection dominated accretion flows, may 
be important in jet formation by creating large-scale 
poloidal magnetic fields [HJ . 
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Therefore one of the most surprising facts related to 
the discovery of NLSyl in 7 -rays is the development 
of a relativistic jet in objects with a relatively small 
BH mass of 10 7 -10 8 M@. However, it is worth noting 
that the mass estimation of the BH in these sources 
may have large uncertainties due to the effect of radi¬ 
ation pressure 36] and the possible disc-like structure 
of their BLR |24j . (]~2l | modelling the optical/UV data 
of some radio-loud NLSyl with a Shakura & Sunyaev 
disc spectrum have estimated higher BH masses than 
those reported in the past, for PMN J0948+0022 and 
PKS 1502+036 comparable to the values estimated for 
blazars. This may solve the problem of the minimum 
BH mass predicted in different scenarios of relativis¬ 
tic jet formation and development, but introduces a 
new problem. How is it possible to have such a large 
BH mass in a class of AGN usually hosted in spiral 
galaxies? Only very sparse observations of the host 
galaxy of radio-loud NLSyl are available up to now. 
Among the NLSyl detected by LAT only for the clos¬ 
est one, 1H 0323+342, the host galaxy was clearly 
detected, suggesting two possible scenarios: the spiral 
arms of the host galaxy [52| or the residual of a galaxy 
merger [MI- Therefore the possibility that the pro¬ 
duction of relativistic jets in these objects could be 
due to strong merger activity, unusual in disc/spiral 
galaxies, cannot be ruled out. 

The accretion rate (thus the mass) and the spin of 
the BH seem to be related to the host galaxy, lead¬ 
ing to the hypothesis that relativistic jets can form 
only in elliptical galaxies [ld|, H3] ■ We noted that the 
BH masses of radio-loud NLSyl are generally larger 
than those in the whole sample of NLSyl (Mbh ~(2 - 
10)xl0 7 M 0 ; [lilHoj), even if still small if compared 
to radio-loud quasars. The larger BH masses of radio- 
loud NLSyl with respect to radio-quiet NLSyl may 
be related to prolonged accretion episodes that can 
spin-up the BHs. In this context, the small fraction of 
radio-loud NLSyl with respect to radio-loud quasars 
could be an indication that not in all the former the 
high-accretion regime lasted long enough to spin-up 
the central BH [47|. 
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